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Abstract 


The  temperature  variations  of  first-,  second-  and  third- soxmd  velocity 

and  attenuation  coefficients  in  one-dimensional  superfluid  helitjm  are 

evaluated  explicitly  for  very  low  temperatures  and  frequencies  (WgT  «  1, 

where  is  the  sound  frequency  and  r  the  characteristic  time) ,  by  using  the 

collisionless  kinetic  equation  and  superfluid  hydrodynamic  equations  with  the 
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energy  dissipation  effect.  The  leading  terms  for  these  sounds  vary  as  T  ,  and 
Che  ratio  of  second  sound  to  first  sound  becomes  unity  as  the  temperature 
decreases  Co  absolute  zero. 
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I.  Introduction 

The  propagation  of  sounds  has  been  an  important  subject  in  the  study  of 

liquid  helixim,  and  especially  the  temperature  variations  of  the  varioiis  sounds 

are  interesting  because  they  are  very  closely  related  to  the  elementary 

excitation  dispersion  relation.  It  is  well  established  that  at  low 

temperatures  and  pressures,  where  roton  excitation  can  be  neglected,  the 

thermal  properties  of  superfluid  helium  are  dominated  by  low*momentum  acoustic 

phonons^  which  do  not  possess  a  normal,  but  rather  an  anomalous  excitation 

spectrum.  These  two  cases  lead  to  different  microscopic  processes.  In  the 

former  case  four-photon  processes  (4PP)  are  dominant,  while  in  the  latter  case 
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3PF  become  Important.  Andreev  and  Khalatnikov  calculated  the  temperature 
variation  of  sound  in  bulk  liquid  helium  using  a  collisionless  kinetic 
equation  together  with  an  equation  for  the  superfluid  velocity.  Singh  and 
Prakash  investigated  the  behavior  of  the  retarded  single -particle  Green's 
function  for  a  weakly- interacting  Bose  gas  and  found  a  new  temperature - 
dependent  term.  However,  they  adopted  the  normal  excitation  spectrum  in  their 
calculations.  Ve  have  obtained  microscopically  anomalous  excitation 
dispersions^  in  two-and  three-dimensional  liquid  helium  which  are  based  on  the 
ring  diagram  approximation,  and  using  these  dispersions^  we  have  not  only 
Improved  on  their  results,^  but  have  also  obtained  first,  second  and  third 
sound  in  thin  helium  films. ^ 

Concerning  second  sound  in  bulk  liquid  helium,  the  velocity  of  this  mode 
is  c^/yj,  where  c^  is  the  velocity  of  first  sound  in  the  low-frequency  limit 

g 

and  low- temperature  phonon  dominant  region.  However,  Maris  has  predicted  a 
new  type  of  second  sound  which  Increases  towards  to  c^  in  the  regime  of 
collinear  phonon  processes.  More  recently,  using  pulsed  time -of- flight  and  cw 

9 

resonance  techniques,  Elsensteln  and  Narayanamurti  have  observed  this  new 
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type  of  collective  mode  as  a  function  of  frequency,  the  so-called  one- 
dlaenslonal  second  sound.  In  superfluid  helium  at  low  temperatures  and 
pressures.  (The  term  "one -dimensional"  means  simply  that  the  thermal  phonons 
In  a  wave  motion  have  wave  vectors  almost  parallel  to  the  propagation 
direction  of  the  wave  In  bulk  superfluid  helliam. )  This  mode  was  predicted  by 
Haris . ^ 

In  many-body  theory  there  are  several  one -dimensional  models^^  of  Bose 
or  Fermi  gases  with  two-body  Interactions  among  particles.  Although  these 
models  can  describe  the  ground-state  energy,  eigenfunction  and  other 
quantities,  they  are  not  appropriate  for  application  to  this  new  mode. 

However,  we  have  recently  developed  a  model^^  of  a  one -dimensional  Bose  liquid 
based  on  the  ring  diagram  approximation.  Through  this  model  we  obtained  the 
pair  distribution  function  of  a  one -dimensional  Bose  liquid  and  then  Its 
anomalous  excitation  spectrum  and  thermodynamic  fimctions.  In  this  mode,  we 
assumed  that  the  Interaction  potential  between  helium  atoms  Is  given  by  a  soft 
potential  with  a  Lennard-Jones  type  tall  as 


^<x)  - 


* 

0 


|x|  <  a 


Ix|  >  a 


(1.1) 


In  this  paper  we  introduce  a  one -dimensional  model  for  superfluid  helium 

to  derive  the  temperature  variations  of  first,  second  and  third  sound  for 
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theory's  sake,  and  also  show  that  the  above  three  sounds  vary  as  T  and  that 
second  sound  reduces  to  first  sound  at  absolute  zero  temperature.  In  the  long 
wavelength  limit,  the  anomalous  excitation  spectrum  is  given  by 
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where  the  coefficients  are 
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where  n  is  nuadier  density.  Ve  note  that  Eq.  (1.2)  represents  an  anomalous 
dispersion  becaxise  B  is  positive.  In  the  case  of  third  sound,  the  excitation 
spectrum  involves  not  only  the  anomalous  excitation  but  also  the  roton 
dispersion,  which  corresponds  to  large  momentum.  For  large  momentum  we  can 
express  the  roton  excitation  spectrum  as 


E(q)  -  ql,^  *  4n.V„ 


(1.4) 


The  energy  E(q)  is  oscillatory.  Around  the  first  minimum  q^,  we  can  replace 
E(q)  approximately  by  a  Landau  form. 


E(q)  -  A  +  (q-q^)' 


(1.5) 
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where  A  and  m*  are  the  energy  gap  and  effective  mass,  respectively. 

In  Sec.  II  we  first  evaluate  the  temperature  variation  of  first  sound 
near  absolute  zero  following  the  method  of  Andreev  and  Khalatnikov,  and  then 
at  low  frequencies  such  that  «  1,  where  is  the  sound  frequency  and  r 
the  characteristic  time,  we  solve  the  superfluid  hydrodynamic  equations  to 
obtain  first  and  second  sound  in  Sec.  III.  In  Sec.  IV  the  derivation  of  one¬ 
dimensional  third  soxind  will  be  given,  and  we  finally  present  numerical 
results  and  discussion  with  a  table  and  graphs  in  Sec.  V.  Throughout  this 
paper  we  take  the  units  such  that  -  1  and  2m  -  1,  where  m  is  the  particle 
mass,  except  for  the  case  in  which  explicit  restoration  of  or  2m  is 
required. 


II.  Kinetic  approach 

The  temperature  variation  of  first  sound  can  be  obtained  by  solving  the 
equations  of  motion  of  the  fluid  assuming  that  roton  excitation  is  negligible 
where  the  phonon  excitation  is  dominant.  In  order  to  obtain  the  sound 
velocity,  we  solve  the  kinetic  equation  for  the  phonon  distribution  function 
n(r,p,t)  without  the  collision  term  together  with  the  equation  of  continuity 
and  the  equation  for  the  superfluid  velocity  v^.  Since  these  equations  were 
given  In  our  previous  paper, ^  we  shall  not  repeat  them  here.  When  the  liquid 
Is  slightly  perturbed  from  equilibrium,  we  can  assume  that  the  deviations  of 
the  distribution  function,  mass  density  and  superfluid  velocity  (represented 
by  n' ,  fi’  and  v'  respectively)  are  small  and  proportional  to  exp[ i(kx-w„t)  ] , 
where  k  is  the  complex  wave  vector,  k^  +  ik2,  which  characterizes  the 
attenuation  of  sound.  Linearizing  the  above-mentioned  three  equations  for  n' , 
p’  and  v^,  wo  obtain 
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where  is  Che  equilibrium  density  of  liquid.  The  nontrivial  solution  can  be 
derived  by  setting  the  corresponding  determinant  equal  to  zero,  and  then 
making  use  of  Eq.  (1.2)  we  obtain  the  sound  velocity  expression 
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where  u  -  (p^/c) (ac/3p^)  is  the  Gruneisen  constant,  and  the  normal  fluid 
density^^  in  Eq.  (2.3)  is  given  by 

B(kgT)^  +  C(kgT)®  +  ...]  .(2.4) 

^0  ^0  ^0 


We  note  that  at  very  low  temperatures  5c  decreases  as  temperature  increases. 


For  low  frequencies  «  1) ,  we  can  apply  the  superfluid  hydrodynamic 

equations  with  the  dissipation  term  to  the  problem  of  first  and  second  sound 
propagation  in  one -dimensional  liquid  helium.  These  equations  are  also  given 
in  our  previous  paper, ^  and  thus  they  will  not  be  discussed  here.  As 
mentioned  earlier,  we  also  assume  that  deviations  of  the  pressure,  temperature 
and  velocity  of  the  normal  component  (F,  T  and  v^)  from  equilibrium  (expressed 
by  P' ,  T'  and  v^)  are  small  and  proportional  to  exp[  i(kx-w^t)  ]  .  Linearizing 
these  hydrodynamic  equations ,  we  obtain  four  hydrodynamic  equations , 
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(3.4) 


where  S  is  the  entropy  per  unit  mass,  and  are  the  mass  densities  of  the 

two  fluid  components,  k  is  the  thermal  conductivity,  and  ^3  ^4 

coefficients  of  second  viscosity.  These  coefficients  must  be  positive,  and 
2 

(Cj^+C2)  5  ^^2^3  satisfied.  Also,  we  must  take  in  view  of 

Onsager's  reciprocal  theorem. 

Evaluating  the  variables  F'  and  T'  from  Eqs.  (3.1) -(3. 2)  and  combining 


them  with  Eqs.  (3.3)-(3.4)  together  with  the  relations 
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we  obtain 
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Equation  (3.7)  is  a  quadratic  equation  In  and  Its  two  solutions 

correspond  to  cwo  nodes  of  sound  propagation  known  as  first  and  second  sound 
with  the  corresponding  attenuation  coefficients.  Solving  Eq.  (3.7)  we  obtain 
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The  expressions  for  first  and  second  sound  and  their  attenuation  coefficients 
are  very  similar  to  those  for  two-dimensional  films. 


IV.  Thir<^  ggun*^ 

When  sound  waves  with  wavelength  longer  than  the  film  thickness 
propagate  along  a  helium  film,  the  normal  fluid  is  clamped  to  the  substrate, 
while  the  superfluid  shows  density  fluctuations.  Since  third  sound  in 

superfluid  helium  films  has  been  observed  in  thick  and  thin  films,  there  are 
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several  theoretical  approaches.  Rutledge,  McMillan,  Mochel  and  Washburn 

(RMMV)  have  analyzed  third-sound  data  in  thin  helium  films  in  terms  of  a 

quantum  hydrodynamic  approach  with  some  phenomenological  assumptions  about  the 

surface  energy  and  roton  excitation.  We  have  also  investigated  thin- film  data 

on  the  basis  of  our  microscopic  theory  without  any  of  the  assximptions  made  by 


RMMW. 


For  the  velocity  of  one -dimensional  third  sound,  we  make  use  of  the 
formula 


03(1)  - 


«(T) 


(4.1) 


where  m  is  the  helium  mass.  Pg(T)  is  the  superfluid  number  density,  and  /c(T) 
is  the  adiabatic  elastic  constant,  which  can  be  obtained  from  the  second 
derivative  of  the  internal  energy  with  respect  to  the  one -dimensional  number 
density.  The  phonon  and  roton  energies  can  be  obtained  from  Eqs.  (1.2)  and 
(1.5)  as 
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and  the  ground  state  energy  is  given  by 
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where  V*  -  (V^  -  6«^/55)  and  7  -  1.3290,  and  for  convenience  we  have  used  the 


following: 
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2  2  2  2  2  2 
Aq  -  D^n  ,  -  D^n  ,  A^  -  D^n 


(4.5) 


Neglecting  the  variation  of  the  roton  energy  A  with  temperature,  we 


arrive  at 
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Thus  we  obtain  third  sound  as 
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where  the  coefficients  are  given  by 
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In  Eq.  (4.7)  the  first,  next  three  and  fifth  terms  are  associated  with  the 
ground  state,  phonon  and  roton  energy,  respectively.  The  remaining  terms  are 
coupled  with  the  phonon  and  roton  energy. 

V.  Numerical  results  and  discussion 

In  this  section  we  present  calculations  based  on  the  formulas  obtained 
in  the  previous  sections.  Before  going  further,  we  first  look  at  the 
excitation  spectrum,  because  all  sounds  we  have  obtained  are  closely  related 
to  the  spectrum.  We  adopted  a  soft  potential  and  obtained  the  excitation 
spectrum  as  a  function  of  the  dimensionless  parameter  qa  shown  in  Fig.  1  with 
the  potential  and  roton  parameters.  As  mentioned  earlier,  our  theory  is 
characterized  by  the  potential  parameters  V^,  and  a.  For  given  values  of 
these  parameters,  we  can  obtain  the  roton  parameters  A  and  q^.  Therefore,  we 
have  determined  the  potential  parameters  such  that  the  excitation  spectnun 
becomes  anomalous  phonon-like  and  roton-like  for  small  and  large  momenta, 
respectively,  and  reproduces  the  same  form  as  that  of  two-dimensional  third 
sound  (Fig.  2)  because  of  the  lack  of  the  experimental  data  in  one -dimensional 
superfluid  helium.  For  the  case  of  helium  films,  we  have  determined  the 
parameters  so  as  to  reproduce  the  specific  heat  data  of  Bretz  et  al,^^  and 
with  this  choice  we  have  obtained  c^  —  84.06  m/s,  which  is  close  to  the 
experimental  value  of  (76  ±  2)  m/s  of  Washburn  et  al.^^  Hence,  the  numerical 
parameters  in  one  dimension  have  been  chosen  in  association  with  the  previous 
values,  which  we  have  taken  earlier  in  two-  and  three-dimensional  liquid 
helium.  They  are  listed  in  Table  1.  However,  the  effective  mass  and  roton 
energy  are  much  lower  and  a  little  higher  than  those  of  chin  helium  films, 
respectively.  The  behavior  of  the  excitation  spectrum,  which  is  analogous  to 
Che  Landau  spectrum,  is  anomalous  phonon- like  for  small  momenta  and  roton-like 
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at  around  for  large  moaenca.  By  increasing  (or  decreasing)  the  numerical 

values  of  n,  a  and  V^,  we  can  make  the  phonon  slope  steeper  (or  slow)  and  also 

deplete  the  roton  minimum  progressively.  By  appropriate  increases  of 

parameters,  the  roton  minimum  may  disappear  and  eventually  become  an 

inflection  point.  Then  the  spectrum  reduces  to  that  of  a  pseudopotential  or 

Coulomic  potential. Neglecting  the  Leonard -Jones -type  tail,  the  excitation 

spectriim  becomes  the  three-dimensional  excitation  spectrum  of  Brueckner  and 
18  19 

Savada.  In  comparison  of  the  excitation  with  Lieb's  results,  ours 
just  corresponds  to  type-1  excitation  (the  so-called  Bogoliubov  spectrum). 

The  thermodynamic  functions  obtained  from  phonon  excitation  (Eq.  (1.2) 

20 

agree  with  Padmore's  results.  Due  to  the  dimensionality,  all  thermodynamic 
functions  are  one  power  lower  in  temperature  than  those  for  two  dimensions. 

For  the  roton  part,  they  have  the  same  form  as  those  for  other  dimensions 
except  for  sisall  modif ications . 

In  Sec.  2  we  have  evaluated  first  sound  (Eq.  (2.3))  by  solving  the 
collisionless  kinetic  equation.  The  leading  term  for  first  sound  increases  as 
T^,  in  contrast  to  T^  and  T^lnT*^  In  two-  and  three-dimensional  liquid  heliian, 
respectively.  In  one  and  two  dimensions  the  absence  of  a  logarithmic  term  is 
due  to  the  fact  that  the  angle  Integrals  become  zero  or  do  not  appear  in  both 
dimensions.  The  first-sound  expression  (Eq.  (2.3))  decreases  monotonically , 
while  first  sounds  obtained  by  considering  of  the  collision  term  as  in  the 
single-collision-time  model^  in  two  and  three  dimensions  increase  gently,  and 
then  decrease  as  temperature  increases  from  absolute  zero.  Ve  note  that  the 
■single-collision- time  model  does  not  hold  in  one  dimension.  The  collisionless 
kinetic  eqxiation  does  not  yield  second  sound,  which  means  that  the  collision 
integral  and  second  sound  are  closely  related  to  each  other. 
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In  the  low-frequency  region,  first-  and  second- sound  propagations  are 
obtained  by  treating  superfluid  hydrodynaoic  equations  including  the  energy 
dissipation  effect.  The  results  are  given  by  Eqs.  (3.9)-(3.12) .  Excluding 
attenuation  terms  and  substituting  thermodynamic  functions  and  fluid 
densities  (see  the  Appendix),  first  and  second  sound  at  absolute  zero 
testperature  reduce  to 

c^<0)  -  CjCO)  -  Aq  -  {2na(V^-|5«^)/m*l‘‘  . 

We  thus  note  that  Aq  in  Eq.  (1.2)  represents  the  sound  velocity  at  absolute 

zero  temperature,  such  that  the  ratio  of  C2  to  c^  becomes  unity  at  absolute 

zero.  We  have  shown  that  the  magnitude  of  second  sound  is  1/72  of  the  first- 

sound  velocity  in  thin  helium  films^  and  find  that  C2  is  1/73  times  that  of 

9 

first  sound  in  the  bulk  case  at  absolute  zero  temperature  in  the  low- 
frequency  region.  Figure  3  illustrates  our  theoretical  sounds  based  on  Eq. 
(3.8)  at  the  density  of  n  -  0.167  A*  ,  which  corresponds  to  2.18  x  10  A  in 
three  dimensions.  The  numerical  value  of  C^(0)  is  given  by  201.4  m/s,  which 
is  close  to  that  of  three  dimensions  rather  than  that  in  two  dimensions.  This 
figure  confirms  that  the  two  sound  modes  show  different  temperature 
variations.  However,  the  overall  variations  of  the  two  sounds  are  very 
similar  to  those  of  two  dimensions.  As  temperature  increases,  second  sound 
passes  through  a  first  gentle  maximum  and  arrives  at  a  minimum  at  around  0.8 
K.  It  then  reaches  another  peak,  after  which  it  decreases  rapdily  and 
vanishes  at  the  critical  point.  This  variation  is  mainly  due  to  the 
thermodynamic  functions  given  in  the  Appendix.  As  temperature  increases  from 
absolute  zero,  first  sound  increases  moderately  and  then  decreaes 
monotonically.  This  indicates  that  the  excitation  spectrum  is  anomalous  at 
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low  monencum .  As  mentioned  in  preceding  sections,  we  note  that  using  a 

phonon-Boltzmann-equation  method,  Maris  shoved  that  second  sound  at  high 

temperatures  and  low  frequencies  increases  gradually  to  first  sound  in  the 

22 

opposite  limit,  and  further  calculations  given  by  Benin  predicted  that  this 

one -dimensional  second  sound  propagates  far  into  the  regime  »  1. 

In  conclusion,  we  have  two  remarks:  (1)  The  leading  terms  in  first, 

2 

second  and  third  sound  vary  as  T  .  (2)  Second  sound  reduces  to  first  sound  at 

absolute  zero  temperature. 
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Appendix 

The  thennod3mamic  functions  and  normal  fluid  density  for  the  phonon 
part  and  roton  density^^  are  given  by 


''ph">  -  k'i  ^  <V>^  ■  *  ■ 


5A^+7 


•  ■  ’  •  ®1  “  2A, 


(IcT)  - 

S.ph^^^  2x‘  A-  .4 


B(kjT)^  +  C(kgT)^  +  . . . ) 


'n.phW  -  U  <V>^  •  '><V>‘  +  <'< V>‘  *  ■■■'  • 


^n.rot^*^^  “  ^2xk^T  ^  %  ® 
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Figure  1.  Elementary  excitation  spectrum  e(q)/kg  for  one -dimensional  liquid 
helium  plotted  against  the  dimensionless  parameter  qa. 

Figure  2.  Temperature  variation  of  third  sound  at  the  density  0.167  A'^. 

Figure  3.  Overall  temperature  variation  of  first-  and  second-sound  velocity 
at  the  density  0.167  A*^. 

Figure  4.  Temperature  dependence  of  the  ratio  density  0.167 

A*^  at  temperatures  below  about  0.3  K. 
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